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Abstract 

We study the spatial dependence of the order parameter near surfaces be- 
longing to the universality class of the ordinary transition. Special attention 
is paid to the influence of a small surface magnetic field h\ at and above 
the bulk critical temperature. A detailed scaling analysis (which is confirmed 
by a perturbative calculation) reveals that h\ may give rise to an anomalous 
short-distance behavior of the order parameter. Close to the surface the mag- 
netization increases with a power law m ~ z K with k = 1 — r]f d ~ 0.23 for 
the three-dimensional Ising model. These results are closely related to ex- 
perimental findings where exponents of the ordinary transition were observed 
in Fe3 Al, while superstructure reflections revealed the existence of long-range 
order near the surface [X. Mailander et al., Phys. Rev. Lett. 64, 2527 (1990)]. 
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A prototypical system to study critical phenomena in restricted geometries is the semi- 
infinite Ising model, terminated by a plane surface and extending infinitely in the direction 
perpendicular to the surface (z-direction) J1J. Spins located in the surface may experience 
interactions different from those in the bulk, for example due to missing neighbors at a free 
surface or due to a strong coupling to an adjacent medium. In the framework of continuum 
field theory such as the 4 model the surface influence is taken into account by additional 
fields like the surface magnetic field h\ and the local temperature perturbation cq at z = 0. 
The latter can be related to the surface enhancement of the spin-spin coupling in lattice 
models [Q]. 

At the bulk critical temperature T c 6 the tendency to order near the surface can be reduced 
(co > 0) or increased (cq < 0), or, as a third possibility, the surface can be critical as 
well. As a result, each bulk universality class in general divides into several distinct surface 
universality classes, called ordinary (c — > oo), extraordinary (c — > -co), and special 
transition (c = c* ). 

Close to the surface, within the range of the bulk correlation length £ ~ l r l _I/ , the 
singular behavior of thermodynamic quantities is markedly changed compared to the bulk. 
For z <C £ the magnetization behaves as ~ jr^ 1 when r = (T — T c b )/T c 6 — > from below, with 
f3\ assuming characteristic values for special and ordinary transition, which are in general 
different from the bulk exponent (3. (At the extraordinary transition the surface is already 
ordered at T c 6 .) Further, the correlation functions near the surface are characteristically 
modified. The correlation function for points within a plane parallel to the surface is given 
by C(Arii) ~ |Ar|||~( rf ~ 2+ ' 7 ll), where the anomalous dimension 7711 is related to (3\ by (3\ = 
{v/2){d — 2 + 7711) [0. Correlations in the z-direction (and all other directions except the 
parallel one) are governed by C(z, z') ~ \z — z'\~( d ~ 2+ri± \ 

Some of the theoretical predictions ]2],|3[] were found in excellent agreement with exper- 
iments carried out by Mailander et al. f4j,||. In these experiments FesAl was studied close 
to the DO3-B2 transition by scattering of evanescent waves generated by total reflection of 
x-rays at a [110] surface. The system was expected to belong to the universality class of 
the ordinary transition, and indeed the exponents measured were in remarkable agreement 
with theoretical predictions ||. A somewhat disturbing feature was that superstructure 
reflections revealed the existence of unexpected long-range order (LRO) near the surface 
reminiscent to the situation at the extraordinary transition. In the sequel it was demon- 
strated by Schmid || that in a similar situation (at the A2-B2 transition in FesAl) an 
effective ordering field hi at the surface can arise when the stoechiometry of the alloy is 
not ideal. Assuming that an h\ is also present at the DO^-B^ transition, the observed LRO 
can be explained, leaving unanswered the question, however, why exponents of the ordinary 
transition were measured despite the LRO. In the following we show that a small hi may 
generate an universal power-law growth of the order parameter near the surface and, as a 
result, a LRO considerably (and, in fact, infinitely) larger than expected from mean-field 
(MF) approximations, while the correlation function is still governed by the exponents of 
the ordinary transition. 

Most of the theoretical studies concerning inhomogenous systems concentrated on the 
behavior at the fixed points Co = ±00 and Co = c* p , respectively. At T c b for both the ordinary 
and the special transition (for h\ = 0) the order-parameter profiles are zero for all z > 0. 
At the extraordinary transition the surface is ordered and the order decays as z~^l v with 
increasing distance from the surface 0, where in the Ising-case /3/u ~ 0.52 ||. Concerning 
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the effects of hi it was assumed for a long time |9j and recently also shown by rigorous 



arguments [10] that a strong hi at the ordinary transition (the so-called normal transition) 
is equivalent to the extraordinary transition. The special transition was studied by Brezin 
and Leibler and by Ciach and Diehl [0 . It was found that at the fixed point the scaling 



field hi gives rise to a length scale l sp ~ h 1 u 1 . For z ^> l sp one finds that m ~ z~^l v as 
at the extraordinary transition. In the opposite limit, z -C l sp , the magnetization behaves 
asm~ mi z^ ~^>l v . Since Pl P < (3, the order still decays, governed by a somewhat smaller 
exponent compared to large distances. For the Ising model (Pl P — (3)/v ~ —0.15 ||. 

What can we expect if a small hi is applied in the presence of a large c , i.e. close to 
the fixed-point of the ordinary transition. In this situation the parameter cq is a so-called 
dangerous irrelevant variable |p|,[13|, comparable to the 4 coupling constant g at and above 
the upper critical dimension d* = 4, and in general must not be naively set to its fixed 
point value cq = oo. Setting the bulk magnetic field h = 0, the remaining linear scaling 
fields at the ordinary transition are r and hi = hi/c^ 0,^,|I4]]. Hence, the behavior of the 



magnetization under rescaling of distances should be described by 

m{z,r, hi) ~ b~ p/v m{zb- x , rb l/u ', hib xi ), (1) 

where the scaling dimension of hi is given by Xi = (d — r]^ rd )/2 0. 

Let us first discuss the profile for h x = 0. As mentioned above, for r > we have m = 
everywhere. For r < 0, on the other hand, the magnetization approaches its bulk value 
nT'b ~ |r| /3//iy for z — > oo. Close to the surface (z £) the magnetization increases with a 
power law [jnj. To see this from (ffl), we set hx = and fix the arbitrary rescaling parameter 
b by setting it ~ z. Then the magnetization takes the scaling form 

m(z,r) ~ z-V" M T (z/£) . (2) 



Since we expect that m(z — > 0) ~ mi and know that mi ~ |r|^° rd , we conclude for the 
short- distance form of the scaling function A4 T ~ Q@° r l v and, in turn, the behavior of m is 
given by m(z,r) ~ \r\^° r z^° r ~P)/ U [|15j. 

We now turn to the case r = and hi ^ 0. This is the situation we are actually interested 
in and which is important for understanding the experimental results of Ref. M. In this 
case the scaling form derived from ([!]) is 

m(zM) ~ z- p/v M hl (zbi /X1 ) . (3) 

First of all we notice from (||) that the scaling field hi gives rise to a length scale l ord ~ h : 1 ^ Xl 
quite comparable to the situation near the special transition discussed above. In order to 
find the short-distance behavior of Ai^i (C) we have to recall that the surface is paramagnetic 



at the ordinary transition. Thus mi will respond linearly to hi ||16|| . Arguing again that 
m(z — >• 0) ~ mi, we now find that M.^ ~ ( Xl for ( < 1, and, in turn, with the scaling 
relation r/± = (7711 + r/)/2 the short- distance behavior is given by 



m(z, hi) ~ hi z K with « = 1 - v± d ■ (4) 



In the opposite limit, z 3> l ord , the magnetization approaches the bulk equilibrium value 

r /pi'n as ~ 7 

-13/u 
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Eq. (|4]) is the central result of this Letter. It states that the magnetization even at (or 
slightly above) T c b in the presence of a surface field hi shows a power-law increase reminiscent 
of the situation below T c 6 . The short-distance exponent k defined in is zero in MF theory. 
Below d*, however, as for the Ising system in d = 3, it is nonzero and positive. Taking the 
literature values for the surface exponents from Ref. 0, one obtains n ~ 0.23, which implies 
a rapid growth of LRO with increasing z. 

The spatial variation of the magnetization discussed above strongly resembles the time 
dependence of the magnetization in relaxational processes at the critical point. If a system 
with nonconserved order parameter (model A) is quenched from a high-temperature ini- 
tial state to the critical point, with a small initial magnetization the order parameter 
behaves as m ~ vn> 1 ' t e , where the short-time exponent 9 is governed by the differ- 



ence between the scaling dimensions of initial and equilibrium magnetization [IB]. Like the 
exponent k in (|4]), the exponent 9 vanishes in MF theory, but becomes positive below d*. 

There is also a heuristic argument for the growth of LRO near the surface. As said above, 
hi generates a surface magnetization mi ~ h\. Regions (on macroscopic scales) close to the 
surface will respond to this magnetization by ordering as well. How strong this influence is 
depends on two factors. First, it is proportional to the correlated area in a plane parallel 
to the surface in a distance z. While in the surface correlations are suppressed, close to 
the surface the correlation length in directions parallel to the surface, £n, grows as ~ z 
|12]| . Second, for small hi (and thus small surface magnetization) it depends linearly on the 
probability that a given spin orientation "survives" in a distance z from the surface. The 
latter is governed by the perpendicular correlation function C(z) ~ z~( d ~ 2+ri ± \ Taking the 
factors together, we obtain 



m[z) 



hxC^tt 1 = h x z x -^\ (5) 



Qualitatively speaking, the surface when carrying a small mi induces a much larger mag- 
netization in the adjacent layers, which are much more susceptible and respond with a 
magnetization m(z) ^> mj. This effect is not observed on the MF level since there the in- 
crease of the correlated surface area is exactly compensated by the decay of the perpendicular 
correlations. 

In order to corroborate our scaling analysis and the heuristic arguments from above, we 
carried out a one-loop calculation for the 4 -model employing the e expansion. Expanded 
in powers of the coupling constant, the magnetization can be written in the form m = 



(0) 



+ gm^ + 0(g 2 ), where mP' is the well known MF solution P,[19"H and m^ is the 
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1-loop term. The latter was calculated exactly for arbitrary cq and hi in Refs. ]TT| , |T2" 
However, the renormalization of uv-divergences and the following improvement with the 
help of the renormalization group where done at (or in the vicinity of) the special transition 
in these references. As a consequence, the anomalous short-distance behavior at the ordinary 
transition was missed. 

The MF solution that satisfies the boundary condition d z m — cm\ z=0 = hi at the surface 
is given by 

m(°>(z) = </Hl ( 6a ) 
V 9 z 

with 



4 



z = z + z + and z , 



cg + 4/n J^/12 



1/2 



Co 



(6b) 



which holds for general Co and ft,].. Close to the ordinary transition (large cq) the mean 
field length scale becomes z + ~ l ord = (12/g) 1 ^ 2 cq/Hx. As expected from (f|), there is no 
anomalous short- distance behavior on the MF level. The profile has its maximum value at 
z — 0, and for z 3> l ord the profile decays as ~ z~^l v with the MF value j3/v — 1. 



The one-loop term rn} 1 ' is given by |12,p0 



m W (z) 



2 Jo 



dz' C(0;z,z')m io \z') I C(p;z',z) 



(7) 



where m^(z) is the zero- loop (MF) profile ( |5a|) and J = (27T) 1 01 J" cf d l p. The propagator 
C(p; z, z') is Fourier-transformed with respect to the spatial coordinates parallel to surface. 



It can be calculated exactly [|lT| , |i2|1 , and the somewhat lenghty results will be omitted here. 
The integrations in (|7|) necessary to obtain the full scaling function A4hi are complicated and 
can only be carried out numerically. However, it is straightforward to extract the divergent 
terms, poles ~ 1/e in dimensional regularization, and the short- distance singularities ~ logz, 
which, when exponentiated, give rise to power laws modified compared to the MF theory. 
Collecting these terms, mS l > is given by (0) with 



C(p; z, z) 



K, 



d-l ~_2H 
Z 



J™ dyy 1 ^ (e 2 ^/ 5 - l) (l + ^ + 



+ finite 



(8) 



where Kd = 2/ (^(AtcY^T (d /2)j and "finite" stands for terms which are finite for e — >• and 
z —>■ 0. Terms of O(l/co) are also omitted in (R|). The zero-momentum propagator C(0, z, z') 



appearing in ((?]) (for cq — > oo) takes the simple form 



C(0,z,z') 



1 1 



5 z\ Zy 
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z 5 



(9) 



where < (>) denotes the smaller (larger) of z and z' . 

Further analysis shows that the uv-divergences can be absorbed in the standard fashion 
by renormalization of the coupling constant K^go = m(1 + 3m/2 + 0(u 2 )) and of the scaling 
field hi )0 = hi(l — w/4e + 0(u 2 )) [0. After this the coupling constant is set to its fixed point 
value u* = 2e/3. Eventually, after exponentiation of logarithms we find the asymptotic 
power laws 



m(z, hi) 



z -l+e/2 for z > yrd 

hi z^ for z < l ord 



(10) 



As expected the decay of the profile for z 3> l ord is governed by the one- loop result (3/u — 
1 — e/2. The short-distance behavior is consistent with our scaling analysis; in first-order e 
expansion k — 1 — i] c £ d = e/6 . 
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Fig. 1: Qualitative shape of the scaling function .Mhi(C) = C /3// ^-^ / 'hi(C) °f the magnetization. 
More details are described in the text. 



A more detailed account concerning the behavior of the magnetization in between the 
asymptotic regimes of Eq. ( [TD[ ) will be given elsewhere [F2I[] . A qualitative preview on the form 
of the scaling function A^hi(C) = C _/3 ^-Mhi(C) ( see Eq. (j|)) is shown in Fig. 1, where the 
asymptotic power laws are quantitatively correct but the crossover is described by a simple 
substitute for the scaling function. Regarding the crossover between ordinary (hi = 0) 
and the extraordinary (or normal) transition (hi = oo) the following scenario should hold. 
While at the ordinary transition m(z) vanishes everywhere, for hi ^ the magnetization 
increases as ~ z K up to z ~ l ord ~ h x 1 ^ Xl (1/xi ~ 1.33 for the Ising model) and thereafter 
crosses over to the long-distance form given in ([T0|) . When hi becomes larger, the short- 
distance increase is steeper and l ord shrinks. For hi — > oo we have l ord — > 0, and one finds 
m ~ z~^l v for all (macroscopic) distances, the result at the extraordinary transition. Largely 
analogous results — monotonous behavior at the fixed points and profiles with one extremum 
in the crossover regime — were found by Mikheev and Fisher for the energy density of the 
two-dimensional Ising model . 



In conclusion, we studied the effects of a small surface magnetic field at the ordinary 
transition at or above T c 6 . We found order-parameter exhibits anomalous short- distance 
behavior in form of a power-law increase m ~ hiz 1 ^^ near the surface implying a much 
more pronounced long-range order than could be expected from mean-field theory. For z 
smaller than l ord and £, the magnetization increases, while the correlation function (and 
thus the structure function) is governed by the asymptotic form of the ordinary transition 
up to correction of O(l/c ) P5| . Assuming that there exists a small ordering field h x in the 
system studied by Mailander et al. ||, our scenario gives a plausible explanation for the 
experimental findings. 
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